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Abstract 

We have simulated structure and dynamics of water in the grooves of a DNA duplex using 
moleculear dynamics simulations. We find signatures of a dynamical transition in both translational 
and orientational dynamics of water molecules in both the major and the minor grooves of a DNA 
duplex. The transition occurs at a slightly higher temperature (Tql ~ 255 K) than the temperature 
(Tl ~ 247 K) where the bulk water is conjectured to undergo a dynamical transition. Groove 
water, however, exhibits markedly different temperature dependence of its properties from the 
bulk. Entropy calculations reveal that the minor groove water is ordered even at room temperature 
and the transition at T ~ 255 K can be characterized as a strong-to-strong dynamical transition. 
The low temperature water is characterized by pronounced tetrahedral order, as manifested in the 
sharp rise near 109° in the 0-0-0 angle distribution. We find that Adams-Gibbs relation between 
configurational entropy and translational diffusion holds quite well when the two quantities are 
plotted together in a master plot for different region of aqueous DNA duplex (bulk, major and 
minor grooves) at different temperatures. The activation energy for the transfer of water molecules 
between different regions of DNA is found to be independent of temperature. 
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I. INTRODUCTION 



Water in the natural world is often found under constrained and/or restricted environ- 
ments, in the hydration layer of proteins and micelles, within reverse micelles and microemul- 
sions, in the grooves of DNA duplex, within biological cells, to name a few. Properties of 
water under such constrained conditions can be quite different from those of bulk, neat wa- 
ter However, it is likely that even under such restricted conditions water retains some of 
its unique properties. Study of these unique properties of water, especially in the hydration 
layer of biomolecules, particularly of proteins and DNA Q, Q, S, H Q, Q , 

has been a subject of great interest in recent times. Hydration layer not only provides the 
stability of the structure of the biomolecules, but also plays a critical role in the dynamic 



control of biological activity. The intercalation of anti-tumor drugs, such as daunomycin, 



into DNA involves active participation of water molecules in the grooves 



13 



The low temperature (near 200 K) "glasslike" transition of hydrated protein has drawn 
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a g reat deal of attention in both experimental and computer simulation studies [15 
[la, [l9l|. Above this transition temperature proteins exhibit diffusive motion and below this 
temperature the proteins are trapped in localized harmonic modes. An important issue in 



recent times is to determine the effects of hydration water on this dynamical transition 20 , 
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24|. Recent studies have shown that dynamics of water in the hydration layer of 



a protein also exhibits strong temperature dependence around the same temperature and it 
seems to undergo a fragile-to-strong transition which preempts an otherwise expected glass 



transition at a lower temperature 
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Study of DNA hydration layer has recently indicated interesting dynamical behavior of 
water in the grooves EJ. Several recent studies have discussed about the origin of the 
slow component of the solvation dynamics in DNA hydration layer [11] . However, a detail 
discussion of this upcoming issue is beyond the scope of this paper. 



A recent computer simulation study by Stanley and co-workers has shown that the liquid- 
liquid (L-L) transition in water induces a dynamic transition in DNA which has striking 
resemblance with that of liquid to glass transition [28]. This study, however, did not explore 
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the dynamics of water in the grooves of DNA. 



There are many questions that have remained unanswered regarding dynamics of groove 
water at low temperature. For example, is there a dynamic transition in the grooves of 
DNA near the L-L transition of bulk water? Does it in any way resemble the one in protein 
hydration layer? Note that the remarkable properties of bulk water have recently been 
attributed to a highly interesting L-L transition at around Tl ~ 247 K, that is, only 26°C 



below the freezing temperature 29|, |3Cj, |3l) . The effects of the bulk water L-L transition on 



groove water dynamics have not yet been investigated. 

In this article we report our finding that water in the grooves of a DNA duplex shows 
a dynamical transition at a temperature (Tgl) slightly higher than the temperature (Tl) 
where the bulk water undergoes the L-L transition. However, the nature and manifestation 
of the transition in the grooves are quite different from that in the bulk. 



II. SYSTEM AND SIMULATION DETAILS 



The system we studied consists of a Dickerson dodecamer DNA duplex (CGC 



GAATTCGCG) |32j solvated in 1565 TIP5P water molecules [23]. We have studied the 
DNA-water system at constant pressure P = 1 atm, at several constant temperatures (NPT 
ensemble) in a simulation box with periodic boundary condition. The molecular dynamics 



simulations of this aqueous DNA system were performed using the AMBER Force Field 



34) 



We have identified the groove water by using the following procedure. We have calculated 
the radial distribution function (g(r)) of water molecules in the system from the major and 
minor groove atoms. On the basis of this g(r), a cut-off distance of 3.5A (the first minima 
of g(r)) from the groove atoms is used for the selection. For bulk water analysis, we have 
considered those water molecules which are beyond 15A from any DNA atoms. We have 
checked that at 15 A away, water indeed regains bulk-like behavior. 
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FIG. 1: Mean square fluctuation (MSF) of DNA duplex (left panel) and diffusivity (right panel) of 
the oxygen atoms of all water in the system. In the left (DNA) panel, MSF of DNA shows a dynamic 
transition at T ~ 247 K. In the right (water) panel, water shows dynamical crossover around same 
temperature from a high temperature power law behaviour (cyan) to a low temperature Arrhenius 
behaviour (red). 

III. RESULTS AND DISCUSSIONS 



A. Mean square fluctuation of DNA and translational diffusivity of water 



We first report the calculated mean square atomic fluctuation (X 2 ) of the DNA atoms 
starting from 300 K down to 210 K in order to characterize the macromolecular "glass" 
transition temperature (Ton a)- Left panel of Figure 1 displays the same. We find that 
the mean square fluctuation (MSF) of DNA slows down dramatically around 247 K and 
continues to remain slow for the lower temperatures. The onset of the change in slope (near 
Tdna ~ 247 K) of MSF indicates a macromolecular dynamic transition. Right panel of 
Figure 1 shows temperature dependence of the diffusivity for all the water molecules in 
the system. It shows a crossover around the same temperature (Tl ~ 247 K) from a high 
temperature power law form to a low temperature Arrhenius form. From the power law fit 



to the high temperature region we get a glass transition temperature o: 
agreement with the earlier simulation study by Stanley and co-workers 



231 K which is in 
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B. Intermediate scattering function 



We next discuss the self intermediate scattering function (ISF) of the oxygen atoms of 
the water molecules in the grooves of DNA for a set of temperatures starting from 300 K to 
210 K. The self intermediate scattering function is defined as 

* _ i , , , . . „ .... / sin\k\ \r(t) — r(0)\\ 
F s (k, t) = (exp(-ik ■ (r(t) - r(0))) = ( |fc||r(t) - r(0)| / ' W 

where k is the wave vector and r(t) is the position of the oxygen atom of the water molecules. 
The \k\ value taken here is 2.5A~ 1 . The translational relaxation time (tt) is obtained by 
fitting the two step relaxation of ISF at different temperatures using Relaxing Cage Model 
(RCM) 



351 ] . The fitting equation used here is given by 

F s (k, t) = [1 - A(k)]e- {t/Ts)2 + A(k)e~ it/TT) ^ (2) 

Here A(k) is Debye- Waller factor, being the translational relaxation time and /3 is the 
stretched exponent. 

Figure 2(a) shows ISF of oxygen atom for the water molecules in bulk, major groove and 
minor groove at 300 K and 260 K. It is evident from both the figures that water molecules 
in both the major and the minor groove tend to behave like a liquid at a temperature lower 
than the bulk. The behavior is more prominent for water molecules in the minor groove. 
This can be ascribed to the fact that translational motion of water molecules in the minor 
groove is more constrained owing to the more ordered structure in the minor groove than 
water molecules in the major groove of DNA. Water molecules in major groove are, in turn, 
translationally more constrained than bulk water. Figure 2(b) shows the temperature 
dependence of for water molecules in bulk, major and minor grooves. For both bulk 
and major groove water the temperature dependence at high temperature region can be 
fitted to the Vogel-Fulcher-Tammann (VFT) law, = t exp[DTo/(T — To)], where D is a 
constant measuring fragility and To is ideal glass transition temperature at which relaxation 
time diverges. In reality, however, the divergence is avoided as below a certain characteristic 
temperature, the functional dependence of relaxation time switches over to an Arrhenius 
form which is a signature of a strong liquid. The crossover temperatures for bulk water 
and major groove water are found to be T L « 247 K and T GL « 255 K, respectively. The 
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FIG. 2: (a) Intermediate scattering function (Fs(k,t)) of the oxygen atoms of the water molecules 
in bulk, major and minor grooves of DNA duplex at two different temperatures, T = 300 K (left 
panel) and T = 260 K (right panel) for \k\ = 2.5A~ 1 .(b) Translational relaxation time (tt) for 
bulk (left panel), major groove (middle panel) and minor groove (right panel) water. Bulk and 
major groove water show dynamical crossover between high temperature VFT behaviour (cyan) 
and low temperature Arrhenius behaviour (red). Minor groove water shows transition between two 
Arrhenius behaviours. 

dynamical transitions are of fragile-to-strong type, although the fragility of major groove 
water is smaller of the two. 

The minor groove water molecules however show a remarkably different translational 
dynamics. Minor groove water does not show any signature of a fragile liquid in the temper- 
ature range studied. Instead, temperature dependence of translational relaxation time for 
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minor groove water fits well into two Arrhenius forms of different slopes with a cross-over 
temperature around Tql ~ 255 K (same as major groove water). This can be understood in 
the context of difference in the structure of hydration layer in the grooves of DNA. Hydra- 
tion in the minor groove is more extensive, regular with a zig-zag spine of first and second 
shell of hydration where as hydration in major groove is restricted to a single layer of water 
molecule 36j . Water molecules in minor groove are thus more structured in comparison with 



major groove water which results in a strong liquid type of behavior for water molecules in 
minor groove even in the high temperature region. This explains why in contrast to bulk 
and major groove water, minor groove water shows a strong-to-strong type of dynamical 
transition. 



C. Orientational dynamics 



We next analyze orientational (dipole-dipole) time correlation function (TCF) of water 
molecules in the different regions of aqueous DNA and the TCF calculated is defined as 

c ^ ~ (MO) -MO)) (3) 

where fi(t) is the dipole moment unit vector of the water molecule at time t and the angular 
bracket corresponds the ensemble averaging. 

Figure 3(a) displays C M (£) for water molecules in bulk, major and minor grooves at 300 
K and 260 K, respectively. Similar to the translational motion, rotation of the minor groove 
water molecules is found to be the most constrained. Figure 3(b) shows the temperature 
dependence of orientational relaxation time (tr) as obtained from stretched exponential fit- 
ting at long time of dipole-dipole TCF for bulk, major and minor groove water. Bulk water 
shows a fragile-to-strong transition around the same temperature (T L ps 247K) as observed 
for translational relaxation time. However, unlike translational relaxation, orientational re- 
laxation shows a transition between two Arrhenius forms of different slopes with a crossover 
temperature Tql ~ 255 K for both major and minor groove water. Strong-to-strong transi- 
tion observed in the minor groove (both translational and orientational dynamics) can be 
attributed to the effect of confinement in the minor groove (higher depth and lower width). 
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FIG. 3: (a) Dipole-dipole time correlation function (C«(i)) of water molecules in bulk, major 
groove and minor groove of DNA duplex at two different temperatures, T = 300 K (left panel) and 
T = 260 K (right panel), (b) Rotational relaxation time (tr) for bulk (left panel), major groove 
(middle panel) and minor groove (right panel) water. Bulk water shows a crossover between high 
temperature VFT behaviour (cyan) and low temperature Arrhenius behaviour (red). Major groove 
and minor groove water show transition between two Arrhenius behaviour. 



23, 



38| and 



It is known that a confined liquid is comparatively less fragile than in the bulk 
this eventually gives rise to a Arrhenius temperature dependence of the relaxation times 
(signature of strong liquids) even at the higher temperature region. The reason for the dif- 
ferent behavior of major groove water is thus probably due to the fact that rotation probes 
local environment more faithfully than translation. 
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IV. MICROSCOPIC CHARACTERIZATION: O-O-O ANGLE DISTRIBUTION 



To understand how the structure of groove water changes across the dynamical transition, 
we have calculated the 0-0-0 angle distribution inside the first coordination shell of a water 
molecule. Angle distributions at three different temperatures (300K, 250K and 230K) for 
groove water molecules are displayed in Figure 4. At all the temperatures, the distribution 
has a two peak character. While the peak at lower angle is the signature of the presence of 
interstitial water molecules inside the first hydration shell, higher angle peak characterizes 
the degree of tetradehrality present. As it is evident from this Figure, with decreasing 
temperature the degree of tetrahedrality increases (higher angle peak height increases) with 
the removal of interstitial water molecules (lower angle peak height decreases) from the first 
hydration shell. Structural change of this kind with decreasing temperature and increasing 
order (as discussed further below) responsible for the dynamical transition for groove water 
molecules. 

Bulk water also exhibits a dynamical transition near 250K. The signatures are, how- 
ever, weaker in the case of bulk water than what are observed in the grooves. The role 
of confinement in fostering the transition of tetraedral water can be understood in the fol- 
lowing fashion. In the confined state, water molecules gain in energy but lose entropy (see 
next section). The tetrahedrally coordinated water is a low entropy state of the system. 
Confinement thus favours the crossover /transition to the tetrahedral state. 



V. ENTROPY CALCULATION 



In order to understand the origin of the large observed differences between the dynamics of 
water molecules in the minor groove and in the bulk, we have calculated the entropy of water 
molecules in the respective regions Q, [3^] at two different temperatures (300 K and 280 K). 
In both the temperatures, minor groove water molecules have substantially lower entropy 
than bulk. At 300 K the difference is ~ 60% of the latent heat of fusion of bulk water. 
Entropy is usually found to be closely correlated with diffusion coefficient and structural 
relaxation time, in Figure 5, we show the correlation between TSconf and translational 
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FIG. 4: 0-0-0 angle distribution of the groove water molecules inside the first coordination shell 
at 300K, 250K and 230K . Note the decrease of lower angle (interstitial) peak and increase of higher 
angle (degree of tetrahedrality) peak height with decreasing temperature. 

diffusivity and show that the Adam-Gibbs relation remains valid for the different regions of 
DNA. Interestingly, we find that the Adams-Gibbs plot for the two different temperatures 
collapse on a single curve which can be fitted to a straight line. This seems to indicate 
that the activation energies for the transfer of water molecules between different regions 
of aqueous DNA are at most weakly dependent on temperature above the L-L transition. 
Dynamics below the L-L transition is too slow to allow a comprehensive study. The present 
calculation of entropy is semi-quantitatively reliable as the entropy of bulk water is correctly 
(within 5 %) reproduced and also the chemical potentials of bulk water and groove water 
are found to be the same, as expected for systems in equilibrium. 

In a thermodynamic co-existence between two phases, a discontinuous change in entropy 
signals the presence of latent heat and a first order phase transition. However, in the present 
case, the large difference in entropy between bulk and minor groove water molecules should 
be regarded as a signature in the difference in structure between the two phases. Because 
of the small number (~65 for major groove and ~30 for minor groove) of water molecules 
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FIG. 5: Adam-Gibb's plot of translational diffusivity (ln(Z>r)) vs 1/TSconf f° r water molecules 
in the different regions of DNA duplex (major groove, minor groove and bulk) at two different 
temperatures, T = 300 K and T = 280 K. 

present in the groove region, a detailed quantification of microscopic structural arrangement 
is hard to perform. 

Because the numbers of water molecules in the two grooves of the dodecamer are rather 
small, we have also simulated a large system with a standardized 38 base pair DNA and 8,000 



water molecules interacting with TIP3P potential [40j. This system is known to sustain a 
stable double helix over a long time period 401 ] . Interestingly, we obtained qualitatively 
similar results for the groove water molecules, but transitions (around 245 K in the grooves) 
are not as prominent since TIP3P is known not to be a good network forming liquid and 
the L-L transition is largely suppressed in the bulk phase. Nevertheless, we do find similar 
kind of transition in the grooves in two different systems with two different water models 
which strengthens the generality of the results obtained in the present study. 



11 



VI. CONCLUSIONS 



In conclusion, we have studied both translational and rotational motions of water in 
the grooves of a DNA duplex . We find that groove water shows a remarkable dynamical 
transition which can explain the transition observed in DNA duplex itself. The fact that 
this transition occurs at not too deeply supercooled water (Tql ~ 255 K) suggests that this 
can be of importance in natural world. 
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